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If gravitation is propagated by a massive field, then the velocity of gravitational waves (gravitons) 



ill depend upon their frequency as (vg/c) 



(c/fXg) , and the effective Newtonian potential 



will have a Yukawa form oc r~^ exp(— r/Ag), where Xg = h/mgC is the graviton Compton wavelength. 
In the case of inspiralhng compact binaries, gravitational waves emitted at low frequency early in 
the inspiral will travel slightly slower than those emitted at high frequency later, resulting in an 
offset in the relative arrival times at a detector. This modifies the phase evolution of the observed 
inspiral gravitational waveform, similar to that caused by post-Newtonian corrections to quadrupole 
phasing. Matched filtering of the waveforms can bound such frequency- dependent variations in 
propagation speed, and thereby bound the graviton mass. The bound depends on the mass of the 
source and on noise characteristics of the detector, but is independent of the distance to the source, 
except for weak cosmological redshift effects. For observations of stellar-mass compact inspiral using 
ground-based interferometers of the LIGO/ VIRGO type, the bound on \g is of the order of 6 x 10^^ 
km, about double that from solar-system tests of Yukawa modifications of Newtonian gravity. For 
observations of super-massive black hole binary inspiral at cosmological distances using the proposed 
laser interferometer space antenna (LISA), the bound can be as large as 6 x 10^® km. This is three 
orders of magnitude weaker than model-dependent bounds from galactic cluster dynamics. 

PACS numbers: 04.80.Cc, 04.30.-w, 97.60.Jd, 97.60.Lf 
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I. INTRODUCTION 

The detection of gravitational radiation by either laser 
interferometers or resonant cryogenic bars will, it is 
widely stated, usher in a new era of gravitational-wave 
astronomy Q|. Furthermore, according to conventional 
wisdom, it will yield new and interesting tests of gen- 
eral relativity (GR) in its radiative regime. These tests 
are generally based on three aspects of gravitational ra- 
diation: its back-reaction on the source, its polarization, 
and its speed. 

(i) Gravitational back-reaction. This plays an im- 
portant role only in the inspiral of compact objects. The 
equations of motion of inspiral include the non-radiative, 
non-linear post-Newtonian corrections of Newtonian mo- 
tion, as well as radiation back-reaction and its non-linear 
post-Newtonian corrections. The evolution of the orbit 
is imprinted on the phasing of the inspiral waveform, 
to which broad-band laser interferometers are especially 
sensitive through the use of matched filtering of the data 
against theoretical templates derived from GR. A num- 
ber of tests of GR using matched filtering of binary in- 
spiral have been discussed, including putting a bound on 
scalar-tensor gravity measuring the non-linear "tail 
term" in gravitational radiation damping ||^, and test- 
ing the GR "no hair" theorems by mapping spacetime 
outside black holes . A concrete test of gravitational 
back-reaction, albeit at the lowest order of approxima- 
tion, has already been provided by the Binary Pulsar 
PSR 1913-)-16, where the tracer of the orbital phase was 
the radio emission from a pulsar rather than matched 



filtering of gravitational waves Q . 

(ii) Polarization of gravitational waves. In GR, 
gravitational waves come in at most two polarization 
states, independently of the source, while in alternative 
theories of gravity, there are as many as six polariza- 
tions Using a sufficiently large number of gravi- 
tational antennas suitably oriented, it is possible to de- 
termine or limit the polarization content of an incident 
wave, and thereby to test theories. For example, should 
an incident wave be shown definitively to have three 
polarizations, the result would be devastating for GR. 
Although some of the details of implementing such po- 
larization observations have been worked out for arrays 
of resonant cylindrical, disk-shaped, and spherical detec- 
tors @J§1, rather little has been done to assess whether 
the ground-based laser-interferometers (LIGO, VIRGO, 
GEO600, TAMA) could perform interesting polarization 
measurements. The results depend sensitively on the rel- 
ative orientation of the detectors' arms, which are now 
cast (literally) in concrete. 

(Hi) Speed of gravitational waves. According to GR, 
in the limit in which the wavelength of gravitational 
waves is small compared to the radius of curvature of the 
background spacetime, the waves propagate along null 
geodesies of the background spacetime, i.e. they have 
the same speed, c, as light. In other theories, the speed 
could differ from c because of coupling of gravitation to 
"background" gravitational fields. For example, in the 
Rosen bimetric theory jlO| with a flat background metric 
T], gravitational waves follow null geodesies of r], while 
light follows null geodesies of g WW 
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Another way in which the speed of gravitational waves 
could differ from one is if gravitation were propagated by 
a massive field (a massive graviton), in which case, Vg 
would be given by, in a local inertial frame, 



? 4 



(1.1) 



where mg and E arc the graviton rest mass and energy, 
respectively. 

The most obvious way to test this is to compare the 
arrival times of a gravitational wave and an electromag- 
netic wave from the same event, e.g. a supernova. For a 
source at a distance D, the resulting value of the differ- 
ence 1 — Vg/c is 



l-!^.5xl0--f^^) /At 
c \ D J \ \s 

where At is the "time difference" , given by 

At = A<a - (1 + Z)Me , 



(1.2) 



(1.3) 



where At^ and Atg are the differences in arrival time 
and emission time, respectively, of the two signals, and 
Z ~ DHq/c is the redshift of the source, with Hq the 
Hubble parameter. In many cases. Ate is unknown, so 
that the best one can do is employ an upper bound on 
Ate based on observation or modelling. The result will 
then be a bound on 1 — Vg/c. 

If the frequency of the gravitational waves is such 
that hf ^ rUgC^, where h is Planck's constant, then 
Vg/c « 1 — i(c/Ag/)^, where Xg = h/nigC is the graviton 
Compton wavelength, and the bound on 1 — Vg/c can be 
converted to a bound on Xg , given by 



Ag > 3 X 10 



12 



km 



D 



1 



lOOHzy/^ / 
200Mpc / J V/At 



1/2 



(1.4) 



The foregoing discussion assumes that the source emits 
both gravitational and electromagnetic radiation in de- 
tectable amounts, and that the relative time of emission 
can be established (by one means or another) to sufficient 
accuracy, or can be shown to be sufficiently small. 

However, there is a situation in which a bound on the 
graviton mass can be set using gravitational radiation 
alone. That is the case of the inspiralling compact bi- 
nary. Because the frequency of the gravitational radi- 
ation sweeps from low frequency at the initial moment 
of observation to higher frequency at the final moment, 
the speed of the gravitons emitted will vary, from lower 
speeds initially to higher speeds (closer to c) at the end. 
This will cause a distortion of the observed phasing of 
the waves and result in a shorter than expected overall 
time Ata of passage of a given number of cycles. Fur- 
thermore, through the technique of matched filtering, the 
parameters of the compact binary can be measured ac- 
curately |l^ , and thereby the emission time At^ can be 



TABLE I. Bounds on Xg from gravitational-wave obser- 
vations of inspiralling compact binaries, using ground-based 
(LIGO/VIRGO) and space-based (LISA) observatories. 
Masses are in Mq. 



mi 7712 


Distance (Mpc) 


Bound on \g (km) 


Ground-based (LIGO/VIRGO) 




1.4 1.4 


300 


4.6 X 10^^ 


1.4 10 


630 


5.4 X 10^^ 


10 10 


1500 


6.0 X 10^^ 


Space-based (LISA) 
lO'^ 10^ 


3000 


6.9 X 10^'^ 


lO'' 10^ 


3000 


5.4 X 10^** 


10^ 10^ 


3000 


2.3 X 10*** 


10* 10* 


3000 


0.7 X 10^^ 



determined accurat ely. Roughly speaking, the "phase in- 
terval" /At in Eq. (g^) can be measured to an accuracy 
1/p, where p is the signal-to-noise ratio. 

Thus we can estimate the bounds on Xg achievable for 
various compact inspiral systems, and for various detec- 
tors. For stellar-mass inspiral (neutron stars or black 
holes) observed by the LIGO/VIRGO class of ground- 
based interferometers, we have D « 200Mpc, / w lOOHz, 
and /At - p-i « 1/10 |ll. The result is Xg > lO^^km. 
For supermassive binary black holes (10* to 10'' M0) ob- 
served by the proposed laser-interferometer space an- 
tenna (LISA), we have D « 3Gpc, / « 10~'^Hz, and 
/At ~ p-i « 1/1000 fill. The result is Xg > lO^'^km. 

We have refined these crude estimates by explicit calcu- 
lations using matched filtering (Table |). We first calcu- 
late the effect of the frequency-dependent massive gravi- 
ton velocity on the observed gravitational-wave phasing. 
We assume that the evolution of the system, driven by 
gravitational back-reaction, is given correctly by gen- 
eral relativity, apart from corrections of fractional order 



{r/XgY, where r is the size of the binary system; these 
corrections can be shown to be negligible for the cases 
of interest. Including GR post-Newtonian (PN) and tail 
terms (1.5PN) in the phasing, and assuming circular or- 
bits and non-spinning bodies, we determine the accuracy 
with which the parameters of the system can be mea- 
sured ("chirp" mass of the system, reduced mass, fiducial 
phase, and fiducial time), and simultaneously find the ac- 
curacy with which the effect of a graviton mass can be 
bounded (effectively, we find an upper bound on X^^). 
We use noise curves appropriate for the advanced LIGO 
detectors, and for the proposed LISA observatory. It is 
interesting to note that, despite the apparent distance 
dependence in Eq. ( |1.4| ), the bound for a given system is 
independent of its distance, because the signal-to-noise 
ratio, which determines the accuracy of /At, is inversely 
proportional to distance. As a result, the bound on Xg 
depends only on the measured masses of the objects and 
on detector characteristics. The only effect of distance is 
a weak Z-dependence arising from cosmological effects. 
The results for the two kinds of detectors, and for var- 
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ious sources are given in Table |. These correspond to 
bounds on the graviton rest mass of order 2.5 x 10"^^ 
eV for ground-based, and 2.5 x 10^^^ eV for space-based 
observations. 

Can bounds be placed on Ag using other observations 
or experiments? If the graviton is massive, then one ex- 
pects that, in the non-radiative near zone of a body like 
the Sun, the gravitational potential will be modified from 
GM/r to the Yukawa form 



V{r) 



GM 



exp(-r/Ag) , 



(1.5) 



Strictly speaking, such a conclusion would require a com- 
plete gravitational theory of a massive graviton, capa- 
ble of making predictions both in the radiative and non- 
radiative regimes, and which otherwise agrees with ob- 
servation. However, as several authors have pointed out 
|p5|-p^, construction of such a theory is a non-trivial 
question. Thus, in the absence of a well-defined theo- 
retical foundation, we shall make the phenomenological 
assumption that, if the graviton is massive in the prop- 
agation of gravitationa l wa ves, the Newtonian potential 
takes the form of Eq. (1.5), with the same value of Ag. 

With this assumption, one can place bounds on Ag 
using solar-system dynamics. Essentially, the orbits of 
the inner planets agree with standard Newtonian gravity 
(including its post-Newtonian GR corrections) to an ac- 
curacy of order 10"^. Since the observed corrections to 
Newtonian gravity in the limit Ag ^ r go as (r/Ag)^ (it 
is the acceleration, not the potential that is important), 
this implies a rough bound Ag > 10"* astronomical units, 
or 10^^ km. Talmadge et al. surveyed solar system 
data in the context of bounding the range and strength of 
a "fifth force" , a Yukawa term added to Newtonian grav- 
ity. The best bound comes from observations that verify 
Kepler's third law for the inner planets: from observa- 
tions of Mars, we find Ag > 2.8 x 10^^ km. Bounds from 



other planets are summarized in Table IV. Apart from 



the Yukawa potential assumption, this bound is solid and 
model independent. 

Thus the bound inferred from gravitational radiation 
observations of stellar mass compact binary inspiral could 
be twice as large as the solar-system bound, while that 
from supermassive binary inspiral as observed by LISA 
could be 2 X 10^ times larger. 

Some have argued for a larger bound on Ag from galac- 
tic and cluster dynamics |2^Jl^,|l^, noting that the ev- 
idence of bound clusters and of clear tidal interactions 
between galaxies argues for a range Ag at least as large 
as a few Megaparsecs (6 x 10^^ km). Indeed this is the 
value quoted by the Particle Data Group |2^. However, 
in view of the uncertainties related to the amount of dark 
matter in the universe, and the absence of a theory that 
can encompass a massive graviton and cosmology, these 
bounds should be viewed with caution. 

The remainder of this paper provides the details under- 
lying these results. In section II, we study the propaga- 
tion of a massive graviton in a cosmological background. 



to find the relation between emission interval and arrival 
interval. In section III, using the standard "restricted 
PN approximation" , in which the gravitational wave- 
form is expressed as an amplitude accurate to the lowest, 
quadrupole approximation, and a phase accurate through 
1.5PN order [0{v/cY] beyond the quadrupole approxi- 
mation, we determine the effect of graviton propagation 
time on the Fourier transform of the waveform, which 
is the central ingredient in matched filtering. In section 
IV, we calculate the Fisher information matrix and de- 
termine the accuracy with which the compact binary's 
parameters can be measured, including a bound on the 
effect of graviton mass. This approach is a reasonable ap- 
proximation to real matched filtering for Gaussian noise 
and large signal-to-noise ratio. We apply the results to 
specific noise curves and binary systems appropriate for 
ground-based (LIGO/VIRGO) and space-based (LISA) 
detectors. Section V discusses bounds on the graviton 
mass using solar-system dynamics. Henceforth, we use 
units in which G = c = 1. 



II. PROPAGATION OF A MASSIVE GRAVITON 

Because some of the detectable compact binaries could 
be at cosmological distances, we study the propagation of 
a massive graviton in a background Friedman-Robertson- 
Walker (FRW) homogeneous and isotropic spacetime. 
We take the line element to have the form E2l 



ds^ = -dt^ + a^{t)[dx^ + Y?{x){de'^ + 



sm 



(2.1) 



where a{t) is the scale factor of the universe and S(x) 
is equal to Xi sin^ or sinhx if the universe is spa- 
tially flat, closed or open, respectively. For a gravi- 
ton moving radially from emitter x = Xe to detector 
X = 0, it is straightforward to show that the compo- 
nent of 4-momentum — constant. Using the fact that 
TOg = —p°'p^gai3 = — a^'^p'^, where E ^ p^ , together 
with p^ / E = dx/dt, we obtain 



dx 
dt 




-1/2 



(2.2) 



where p^ — a'^{te){E'^ — mg). Assuming that Ee'> mg, 
expanding Eq. ( ^.2| ) to first order in {nig/ Ee)"^ , and in- 
tegrating, we obtain 



Xe 



1 



dt 



i{t) 2a^{te)ElJt 



a(t)dt . 



(2.3) 



Consider gravitons emitted at two different times and 
tg, with energies E^. and E'^, and received at correspond- 
ing arrival times (xe is the same for both) . Assuming that 
A^e = ig — t'e ^ a/d, and noting that vigj E^ = {Xgfe)~^, 
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where /e is the emitted frequency, we obtain, after eHm- 
inating Xe, 



Ata = il + Z) 



D 
2Af 



1 



1 



(2.4) 



where Z = qq / a(te) — 1 is the cosmological redshift, and 



ao 



a{t)dt , 



(2.5) 



where ao = a{ta) is the present value of the scale factor. 
Note that D is not a conventional cosmological distance 
measure, like the luminosity distance = aoS(xe)(l + 
Z), or the proper distance Dp = aoXe- For Z <C 1, it is 
given by the standard formula D = Z/Hq] for a matter 
dominated, spatially flat universe, D and are given 
by 



D = (2/5i/o)(l + Z){1 - (1 + Z)-5/2) _ 
= (2/i7o)(l + Z)(l-(l + ^)-'/'). 



(2.6a) 
(2.6b) 



The ratio D /Dl will play a role in our analysis of the 
bound on Xg. It has the following representative behav- 



D _\ l-Z + 0{Z^), 

I 5{i+zr 



z < 1, aU no 
no = 1, aU z 



(2.7) 



where ilo is the density parameter. At Z = 1, the factor 
D/Dl varies from 0.5 for = 0.01 to 0.6 for flo = 2. 
For simplicity, we shall henceforth assume that JIq = 1- 



III. MASSIVE GRAVITON PROPAGATION AND 
THE PHASING OF GRAVITATIONAL WAVES 

We shall treat the problem of a binary system of com- 
pact bodies of locally measured masses mi and m2 in 
a quasi-circular orbit, that is an orbit which is circular 
apart from an adiabatic inspiral induced by gravitational 
radiation reaction within GR. We ignore tidal interac- 
tions and spin effects. For matched filtering of gravita- 
tional waves using LIGO/ VIRGO or LISA type detec- 
tors, it is sufficient for our purpose to write the gravita- 
tional waveform h{t) in the "restricted post-Newtonian 



form" p3y24yi2[| , in terms of an amplitude A{t) expressed 
to the lowest, quadrupole approximation, and a phase 
^(t), expressed as a post- Newtonian expansion several 
orders beyond the quadrupole approximation. 



h{t) = A(t)e-**(*) , 

$(t) = <i>c + 2nj^ f{t)dt, 



(3.1a) 
(3.1b) 



where f{t) is the observed frequency of the waves, and 
$c and tc are "fiducial" phase and time respectively. The 
amplitude A is given by 



A{t) 



2fj, 



aoS(xe) r(t) 



F(j,0,0,V), 



(3.2) 



where m = mi + m2 and fi = mim2/m are the total and 
reduced mass of the system (we also define the reduced 
mass parameter 77 = /i/m), r(t) is the orbital separation, 
and F is an angular function related to the orientation 
of the orbit (angles i, tp) and the direction of the source 
relative to the antenna (angles 9, 0), given by 



F^{i,0,<j),i^) = ^{1 + cos^ if Fl + cos^ iF, 



2 

X I 



(3.3) 



where F+{9, 0, ip) and Fx {6, (j), -0) are beampattern fac- 
tors quoted, for example in Eqs. (104) of Q|. For sim- 
plicity, we shall average over all four angles, and use the 
fact that (F^) =4/25. 

We next compute the Fourier transform of h{t). Ex- 
panding h{t) about the time t at which the observed fre- 
quency is /, i.e. f{t) = /, and using the stationary-phase 
approximation, we obtain 



where 



4 Me 



5 aoS(xe) 



(^Xe/e)'/' 



(3.4) 



(3.5a) 



^if) = 27r / {t- Qdf + 2nft, - $e " 7r/4 , (3.5b) 

where Ad^ = rf/^m is the "chirp" mass of the emitter, 
and where we have used the Newtonian relation m/r(t) = 
{-Kuifef^^. The subscript "e" denot es "at the emitter". 
We next substitute Eq. (2^) into ( 3.5b ) to relate the 
time at the detector to that at the emitter, noting that, 
because of the cosmological redshift, /e = (1 + Z)/. The 
result is 



«-(/) = 2tI j {te - tec)dfe 

fee 



ttD 



2nfte 



' 4 ' 
(3.6) 

where = - D/[2{1 + Z)\lf% and = 
27r£'/[(l -I- Z)Ag/c]. To find te — tec as a function of /e, 
we integrate the equation for radiation reaction between 
tec and tp-. 



4/e 
dte 



96 



1 



743 11 ^ , ^ ,2/3 
336 4 " ^ 



-47r(7rm/e 



(3.7) 



where we have included the first post-Newtonian (FN) 
term and the 1.5PN "tail" term in the radiation- reaction 
equation (see, e.g. p3|). After absorbing further con- 
stants of integration into tc and dropping the bars 
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on those two quantities, and re-expressing everything in 
terms of the measured frequency / [note that (/)^^^ = 
(d/e/die)^/V(l + Z)], we obtain 



f i(/)e''*(/) 

lo, 




(3.8a) 
(3.8b) 



*(/) = 



27r/tc - $c - 7r/4 + 1*-^/^ 

J c c / ^ 



5 /743 11 N 2/-5 

H — — I — V ir ' 1 

96 V336 4 " ' 



(3.8c) 



where u = nMf, and is the "measured chirp mass", 
related to the source chirp mass by a redshift: M = 
(1 + Z)A4e- The parameter f3 is given by 



AfO + Z) 



(3.9) 



The frequency /max represents an upper cut-off fr equ ency 
where the PN approximation fails. Equations (3^) are 
the basis for an analysis of parameter estimation using 
matched filtering. 

Before turning to matched filtering, we must address 
our approximation of the motion and gravitational radi- 
ation damping as being general relativistic up to correc- 
tions of order (r/Ag)^. In the radiation- reaction formula 



Eq. (3.7), we included corrections to the quadrupole 
formula at 1.5PN order, corresponding to corrections 



of order 



Thus our neglect of massive graviton ef- 



^Xg'^y-^ < 1 for all 
m/r for circular or- 



fects amounts to assuming that 
systems of interest. Because v"^ 

bits, we can rewrite this condition as (m/Ag)u~^/^ <C 1. 
Since typically v > 10~^ for all systems of interest, and 
Xg > 10^^ km from solar-system bounds, this condition 
is easily satisfied. 



IV. BOUNDS ON THE GRAVITON MASS USING 
MATCHED FILTERING 

A. Matched-filter analysis 

To obtain a more reliable estimate of the bound that 
can be placed on the graviton mass, we carry out a full 
matched-filter analysis following the method outlined for 
compact binary inspiral by Cutler and Flanagan |23| and 
Finn and Chernoff [|^. The details here parallel those 

of g. 

With a given noise spectrum Sn{f), one defines the 
inner product of signals hi and by 



(/ll|/l2) = 2 



Snif) 



(4.1) 



where ha is the Fourier transform of the waveform de- 
fined in Eqs. (3^) (henceforth, we drop the tilde on 



frequencies). The signal-to-noise ratio for a given signal 
h is given by 



p[h] = s/N[h] = {h\hy/^ . 



(4.2) 



If the signal depends on a set of parameters 0° which are 
to be estimated by the matched filter, then the rms error 
in 9°' in the limit of large S/N is given by 



(4.3) 



where is the corresponding component of the inverse 
of the covariance matrix or Fisher information matrix Tab 
defined by 



Tab 



/ dh 


dh\ 


I 961° 


36^ I 



(4.4) 



The correlation coefficient between two parameters 
and 9^ is 



(4.5) 



We estimate the following six parameters, In^, $c, 
/o^c, InAl, In 77, and /3, where /o is a frequency charac- 
teristic of the detector, typically a "knee" frequency, or 
a frequency at which Snif) is a minimum. The corre- 
sponding partial derivatives of h{f) are 



dhjf) 
d\nA 
dhjf) 

9$c 

dhjf) 
dfotc 

dhjf) 
dlnM 



dhjf) 
dluT] 



dhjf) 
df3 



Kf), 

2m{f/fo)hif ) , 



128 



96 



7(77)^" 



-^^-3/5^-2/3 



5i 
96 

lo" 

-1 



ri-y'u-y']hif).. 



-lu-'hif), 



(4.6a) 
(4.6b) 
(4.6c) 

(4.6d) 

(4.6e) 
(4.6f) 



where 7(77) = (743/336 llr;/4)r;-2/5^ and 7' = ^7/^7;. 
Since we plan to derive the error in estimating (3 about 
the nominal or a priori GR value /3 = 0, we have set 
/3 = in all the partial derivatives. 

We assume that the detector noise curve can be ap- 
proximated by an amplitude So, which sets the overall 
scale of the noise, and a function of the ratio ///o = 
X, which may include additional parameters, that is 
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TABLE II. The rms errors for signal parameters, the corresponding bound on Ag, in units of lO^'^ km, and t he co rrelation 
coefficients CAirj, cmi3 and The noise spectrum is that of the advanced LIGO system, given by Eq. (4.10), and a 

signal-to- noise ratio of 10 is assumed. Masses are in units of A/q, Ate is in msec. 



mi 


m2 




At, 


AM/M 


Ari/ri 


\ 




Cm 13 




1.4 


1.4 


4.09 


1.13 


0.034% 


7.88% 


4.6 


-0.971 


-0.993 


0.992 


1.4 


10.0 


6.24 


2.04 


0.191% 


12.2% 


5.4 


-0.978 


-0.994 


0.994 


10.0 


10.0 


9.26 


3.53 


1.42 % 


57.3% 


6.0 


-0.983 


-0.994 


0.997 



Snif) — Soga{x), where the sub scri pt a d eno tes a set 
of parameters. Then from Eqs. (3.8) and (4.2) we find 



that the signal-to-noise ratio is given by 
p = 2^/-2/^(/(7)/5o)i/2 



15 Dl 

where we define the integrals 



-q/3 



9a{x) 



-dx . 



(4.7) 



(4.8) 



Note that any frequency cut-offs are to be incorporated 
appropriately into the endpoints of the integrals I{q). 
If we define the coefficients Iq = I{q)/I{7), then all 
elements of the covariance matrix turn out to be pro- 
portional to p^ times linear combinations of terms of 
the form Uq "^'^/g for various integers n and q, where 
uq = nMfo. This overall p dependence is character- 
istic of the large S/N limit. As a result, the rms errors 
Afl" are inversely proportional to p, while the correlation 
coefficients are independent of p. Defining A/? = A^/^/p, 
viewing A/3 as an up per b ound on_/3, and combining this 
definition with Eqs. (3^) and (4/7) we obtain the lower 
bound on A^: 



Xg > 



_2_/(7) 
15 S-o 



1/4 



D 



{1 + Z)Dl 



^2/3_V(ll/ 



12 



(4.9) 



Note that the bound on Ag depends only weakly on 
distance, via the Z dependence of the factor [D/{1 + 
Z)Dl]^^^, which varies from unity at Z = to 0.45 at 
Z = 1.5. This is because, while the signal strength, and 
hence the accuracy, falls with distance, the size of the 
arrival-time effect increases with distance. Otherwise, 
the bound depends only on the chirp mass and on detec- 
tor noise characteristics. We now apply this formalism 
to specific detectors. 



B. Ground-based detectors of the LIGO/VIRGO 
type 

The proposed advanced version of LIGO is expected to 
detect compact binary inspiral to distances of 200 Mpc 



to 1 Gpc. The sensitive frequency band extends from 
around 10 Hz to several hundreds of Hz. We adopt the 
benchmark advanced LIGO noise curve, given by 



Snif) 



00 , 

Soiifo/f)' 



2(///o)']/5. 



/ < lOHz 
/ > lOHz 

(4.10) 



where 5*0 = 3 x 10-''*Hz"\ and /o ^ 70 Hz. The cut- 
off at 10 Hz corresponds to seismic noise, while the 
and dependences denote thermal and photon shot- 
noise, respectively |13 . We choose an upper cut-off fre- 



quency, where the PN approximation fails, correspond- 
ing to the innermost stable circular orbit. Although this 
is known rigorously only for test body motion around 
black holes [Em, a conventional estimate is given by 
/isco ~ [6'^/^(rni + m2)]~^. Converting this to the 
measured frequency and chirp mass, we have Xmax = 
[6^/^7r77^^/^Al/o]"^. For this case, we thus have g{x) — 
(x-"^ + 2 + 2x^)/5, and I{q) = ™[a:-«/Vff(a;)]dx. We 
then calculate and invert the covariance matrix and eval- 
uate the errors in the five relevant parameters (the pa- 
rameter In A decouples from the rest and is relevant only 
for the calculation of p), and the correlation coefficients 
between A4, rj and (3. For various "canonical" compact 
binary systems observable by advanced LIGO, the re- 
sults are shown in Table ||. Note that, in determining 
the bound on Xg, we must include the Z dependence em- 
bodied in Eq. (|4.9| ). To do so, we take our assumed 
value for signal-to-noise ratio p — 10, determine the 
luminosity distance using Eq. (4.7), and convert that 
to a redshift using Eq. (2.6b), with an assumed value 



i/o — 50 km s Mpc and VLq — 1. W e then substitute 

it, along with Eq. (|2.6aD into Eq. (|0). 



It is useful to compare these results to those from pa- 
rameter estimation calculations using pure GR to 1.5PN 
order mcZwdm^ spin-orbit effects (see e.g. [|3 25|). There, 
an additional parameter related to the spin-orbit effect 
(also called /3, with a nominal value of zero) was es- 
timated, although it produced a different u-dependent 
term in the phasing formula (it~^/^ instead of u~^). Nev- 
ertheless, the errors in the fiducial phase A<I>c, time /S.tc 
and chirp mass A In are virtually identical in both 
cases, and somewhat larger than if no additional (3 pa- 
rameter were estimated (compare Table ^ with Table 
I and II of m or Table II of [||). But in our case, 
the errors in the reduced mass parameter 77 are signif- 
icantly larger, a result of the nearly perfect correlation 
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TABLE III. Rms errors on signal parameters, the bound on \g, in units of 10^^ km, and the c orrelation coefficients. The 



noise spectrum is that of LISA including white-dwarf binary confusion noise, given by Eq. (4.11). Signal-to-noise ratio p is 
shown, corresponding to a luminosity distance of about 3 Gpc. Masses are in units of Mq, Ate is in sec. 



mi 


m2 


p 


A(j)c 


Ate 


AM/M 


Arj/r] 


K 

y 








10' 


10'^ 


1600 


0.073 


20.0 


0.0187% 


0.562% 


6.9 


-0.979 


-0.992 


0.997 


10^ 


10« 


710 


0.145 


22.5 


0.0119% 


0.362% 


3.9 


-0.984 


-0.995 


0.997 


lO'' 


10** 


5800 


0.017 


0.48 


0.0021% 


0.108% 


5.4 


-0.954 


-0.985 


0.991 


10'' 


10^ 


4300 


0.026 


0.40 


0.0015% 


0.062% 


3.0 


-0.970 


-0.992 


0.991 


10^ 


10^ 


2100 


0.017 


0.09 


0.0008% 


0.072% 


2.3 


-0.946 


-0.975 


0.992 


10^ 


10* 


750 


0.048 


0.18 


0.0007% 


0.059% 


1.2 


-0.955 


-0.987 


0.989 


10^ 


10" 


320 


0.092 


0.22 


0.0004% 


0.141% 


0.7 


-0.963 


-0.992 


0.989 



{u~^ dependence) b etwee n the IPN term and the /3-term 
in the phasing, Eq. (3.8c). The error grows dramatically 
with total mass because the smaller number of observed 
gravitational-wave cycles reduces the ability of the tail 
term (oc u~^/^) to break the degeneracy. 



For supermassive binaries ranging from 10* to IO^Mq, 
and for integration time T — 1 year, we estimate the 
errors in the five parameters and determine a bound on 
Xg. We choose the signal-to-noise ratio p for each case 
such that the luminosity distance to the source ~ 3 Gpc, 
so that cosmological effects do not become too severe. 
The results are shown in Table III. 



C. Space-based detectors of the LISA type 

The proposed Laser Interferometer Space Antenna 
(LISA) is expected to be able to detect the inspiral of su- 
permassive black hole binaries to cosmological distances, 
with very large signal-to-noise ratio. The sensitive fre- 
quency band extends from around 10^* to 10^^ Hz, with 
a typical integration time of the order of one year. We 
adopt a noise curve described in the LISA pre-Phase A 
report |Q , augmented by a fit to "confusion noise" gen- 
erated by a population of close white dwarf binaries in 
our galaxy |^^, given by the equations: 

5*0 = 4.2 X lO^-^^Hz"^ , 
/o = lO-^Hz , 



(4.11) 



g^x) = VlOx""/^ -I- 1 + xVlOOO 

_|_323,52;-(6-398-l-3.5181ogioa:) 



In order, the four terms in g(x) correspond to: tempera- 
ture fluctuations, photon shot noise, loss of sensitivity 
when the arm lengths exceed the gravitational wave- 
length, and a fit to the white-dwarf binary confusion 
noise. For the maximum frequency, we again adopt that 
of the innermost stable circular orbit. The minimum fre- 
quency is set by the characteristic integration time for 
LISA, nominally chosen to be one year. We calculate the 
time Te remaining until the system reaches the innermost 
stable orbit by integrating Eq. (3.7) using only the dom- 
inant, Newtonian contribution, convert fr om t ime at the 
emitter to observation time T using Eq. (2.4) (ignoring 
the small correction due to massive graviton propaga- 
tion), and obtain 



1 



3/8 



nMfo \256T J 



(4.12) 



V. SOLAR-SYSTEM BOUNDS ON THE 
GRAVITON MASS 

If the Newtonian gravitational potential is modified by 
a massive graviton to have the Yukawa form of Eq. (|l.E|), 
then the acceleration of a test body takes the form 



(5.1) 



where 



lj{r) = M (1 + r/Xg) cxp{-r/Xg) 



M 



(5.2) 



For a planet with semi-major axis Qp and period Pp, 
Kepler's third law gives ap(27r/Pp)^/^ — fi{ap)^^^. For 
a pure inverse-square law, /j, = constant, and its value 
is determined accurately using the orbit of the Earth. 
Thus, by checking Kepler's third law for other planets, 
one can test the constancy of fi. For a given planet, we 
define the parameter -qp by 



1/3 



(5.3) 



Combining Eq. (5_^) and (^.2|), we obtain a bound on Xg 
in terms of rjp 



Xg > 



1 - a, 
6r]p 



2\ 1/2 



(5.4) 



where Xg and Op are expressed in astronomical units 
(1.5 X 10^ km). Table |^ lists the observed bounds on 
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TABLE IV. Bounds on \g in units of 10^^ km from Kepler's 
third law applied to the solar system. Semi-major axes are in 
astronomical units, and the appropriate one-sided, 2a bound 
on rip from Talmadge et al. IISI is shown. 



Planet 


Bp 


Bound on rjp 




Mercury 


0.387 


1.4 X 10"* 


0.5 


Venus 


0.723 


1.5 X 10"^ 


1.1 


Mars 


1.523 


-6.5 X 10"^° 


2.8 


Jupiter 


5.203 


-6 X 10"** 


1.3 



rjp for Mercury, Venus, Mars and Jupiter compiled by 
Talmadge et al. ]l9| , and the resulting bounds on \g. 

A Yukawa violation of the inverse square law will also 
produce anomalous perhelion shifts of orbits, of the form 
5uj ~ TT^Gp/ Xg)'^ . By comparing measured shifts with 
the general relativistic prediction for Mercury and Mars 
(Venus and Jupiter are unsuitable for this purpose) [p^ , 
one obtains bounds a factor two weaker than those from 
Kepler's law. 

Systems like the binary pulsar do not yield useful 
bounds, mainly because they are so compact. Since the 
deviations from GR go as (a/Ag)^, where a is the size 
of the system, the bound is roughly Ag > a/e^/^, where 
e is the accuracy of agreement with observations. Since 
a ~ 10® km, and e ~ 10~^, the bound is far from com- 
petitive with that from the solar system. 
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